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Abstract. A problem to account for the direct electron-positron annihilation in 
positron-hydrogen scattering above the positronium formation threshold has been 
resolved within the time independent formalism. The generalization of the optical 
theorem is derived for the case when an absorption potential is present in the 
Hamiltonian. With this theorem the annihilation cross section is fully determined 
by scattering amplitudes. This allows us to separate out analytically the contribution 
of the positronium formation from the overall annihilation cross section. The rest is 
determined as the direct annihilation cross section. It is done uniformly below as well 
as above the positronium formation threshold. The multichannel three-body theory for 
scattering states in the presence of an imaginary absorption potential is developed in 
order to compute the direct e + e~ annihilation amplitude. Special attention has been 
paid to an accurate definition of the coordinate part of the absorption potential as 
the properly constructed zero-range potential, which corresponds to the delta- function 
originated from the first order perturbation theory. The calculated direct annihilation 
cross section below the positronium formation threshold is in good agreement with 
results of other authors. The direct annihilation cross section computed with the 
formalism of the paper shows nonsingular behavior at the positronium formation 
threshold and is in good agreement with existing data. A number of e + e~ direct 
annihilation cross sections and positronium formation cross sections in the energy 
gap between Ps(ls) and H(n = 2) thresholds are reported. A sharp increase in the 
calculated direct annihilation cross section across the resonant energy is found for all 
first S and P-wave Feshbach resonances. 
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1. Introduction 

In positron-hydrogen collision, if the energy is higher than the positronium formation 
threshold, we have two genuine different asymptotic channels. The positron may go 
away leaving the hydrogen atom behind in ground or in an excited state, or the positron 
and the electron may form a positronium leaving the proton behind. This picture 
gets modified if we take into account the possibility that the positron and the electron 
annihilate. In case of rearrangement scattering, when the positron and the electron form 
positronium in the Ps(ls) state, the positron-electron pair annihilates mostly from this 
state. The lifetime of the Ps(ls) state depends on the total spin pQ. In this process, once 
positronium is formed, the third particle does not affect the annihilation. The direct 
annihilation occurs without formation of the positronium. In this case the positron 
interacts with the electron in the field of the proton only when all the three particles 
are close to each other. Thus, the direct annihilation is a genuine three-body process. 

The standard approach treats the e + e~ annihilation on the basis of the perturbation 
theory [21 [3]. The annihilation cross section appears as a matrix element of 
the coalescence operator between scattering states, which are the solution of the 
unperturbed three-body Coulomb Hamiltonian. At higher energies, this formalism have 
difficulties. At the positronium formation threshold, the calculated annihilation cross 
section becomes infinite. There has been a number of attempts to improve theory and to 
achieve a unified treatment of the positron annihilation and the positronium formation 
[H EJ Ej. In one approach, to make the lifetime of the e + e~ pair finite, an imaginary 
absorption potential is added to the positron-electron subsystem Hamiltonian [I] . " The 
loss of particles", due to this potential, is then interpreted as the e + e~ annihilation. 
The standard QED formula [21 [3] for the positron annihilation cross section suggests 
the coordinate part of the absorption potential to be the three-dimensional Dirac delta 
function depending on the relative position vector between the positron and the electron. 
Subsequently, this kind of absorption potential was used for computing the annihilation 
cross section above the positronium formation threshold [El [7J E] . The three-dimensional 
delta-function is too singular to be an ordinary potential in the Schrodinger equation 
[H [101 CLE]. In the actual calculations in Refs. [H El [71 [8] with the delta potential, this 
singularity was smoothed out. Nevertheless, the mathematically sound formulation of 
using the three-dimensional zero-range potentials in three-body Coulombic systems is 
still needed. It is done in this paper to complete the theory. 

Direct annihilation cross section above the positronium formation threshold is 
subject of an extensive study [12J. In Ref. [HI [6] it was claimed that the direct process 
and the annihilation after the positronium formation are inseparable and therefore the 
notion of the direct annihilation above the threshold is meaningless. In contrast, the 
direct cross section was defined and computed successfully below and above the Ps- 
formation threshold in Ref. [5]. This calculation showed no sign of the presence of any 
threshold behavior in the energy dependence of the direct annihilation cross section. 

This paper is devoted to continue the study of the direct annihilation within the 
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time independent formalism introduced in [IB] . In section 2 we define the scattering 
solution for the e + — H Schrddinger equation in the presence of an absorption potential. 
Subsequently, the multichannel formalism, which is needed to determine scattering and 
absorption amplitudes above the positronium formation threshold, is developed. In 
section 3 we derive the generalization of the optical theorem in the presence of an 
absorption potential. This optical theorem is used to separate out the annihilation after 
the positronium formation cross section from the overall absorption cross section. The 
remaining part is naturally the direct annihilation cross section. The correct form of 
the zero-range absorption potential is introduced in section 3. Section 4 contains the 
results of calculations for the direct annihilation. Section 5 concludes the paper. 

Throughout the paper we use bold letters for vectors, e.g. x, and non-bold for 
their magnitudes, e.g. x = |x|. The unit vector pointing out in the direction of the 
vector y is denoted as y = y/y. In the paper we work with wave-functions and their 
components of the form \l/(X,p) where X stands for the configuration coordinate and 
p for the momentum of the asymptotic motion. For the sake of clarity of notations and 
to avoid overloading of formulae we drop very often the configuration space coordinate 
from the notations of wave-functions and their components leaving only the momentum 
explicitly. 



2. Three-body scattering formalism 

In this section we apply the three-body multichannel formalism to e + — H scattering. 



2.1. Three-body scattering in the presence of absorption potential 

We consider the three-body problem with an additional imaginary potential which acts 
between the positron and the electron. If the center of mass motion is removed, the 
Hamiltonian in the Jacobi coordinate system takes the form 

H = H c + igW 2 (x 2 ), (1) 
H C = H + Vf( Xl ) + Kf (x 2 ) + V 3 c (x 3 ), (2) 



— V 2 - — V 2 V- C (x.-) = e -^- k 



Ho = -^-V 2 - £-V 2 Vf( Xi ) = ^. (3) 



Here we assign numbers 1,2 and 3 to the positron, the proton and the electron, 
respectively, therefore the electric charges are expressed through the unit charge e as 
e\ = e 2 = e, e^ = —e. The Jacobi coordinates are defined in terms of particle radius- 
vectors Ti by the standard formulae 

Xj = Tj - r fe , = ri ■ , (4) 

rrij + m k 

and the reduced masses are given by 

m^mu mArrin + mi) . , 

im = — — —, n = — ■ (5) 

rrij + nik rrii + rrij + m,k 
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The potential igW 2 (^ 2 ) acting between particles of the pair 2 is defined such that g is 
real and negative, and W 2 (k 2 ) is real and non negative. In this case igW 2 is a complex 
absorbing potential. We do not specify the coordinate dependence of W 2 yet, except of 
requiring that W 2 is short-ranged. 

The Schrodinger equation for the positron-hydrogen scattering reads 

(H + Vf - E)$ + = -{V 2 C + Vf + igW 2 )* + . (6) 

The scattering solution is defined at real energy E by the asymptotics as y\ — > oo 

oc MxiW* 1 ' 71 + HPiYuPi)}- (7) 

Here 0i(xi) is the hydrogen ground state wave function with the energy e\, and 
the incident momentum of the positron pi is related to the energy by the condition 
E = h^2vxp\ + e\. We note, that due to the asymptotics of the wave-function (j7J), 
the term V 2 C + ^3 m the right hand side of Eq. ([6]) is always confined into the region 
of the configuration space where the hydrogen wave function 0i(xi) is not negligible. 
This makes the term V 2 C + Vg short-range-type and hence the asymptotics ([7]) in the 
coordinate yi is free from the Coulomb contribution. This property holds true for all 
equations we deal with in the paper. 

Conventionally |14j . the scattering amplitude F can be represented through the 
wave function. It is done by rewriting Eq. ([6]) in the integral form and taking the 
asymptotics y\ — > 00. By doing so, we get the Lippmann-Schwinger equation (LSE) 

tf+( P i) = $x( Pl ) + G 1 (E + )(V 2 C + V 3 C + igW 2 )^( Pl ), (8) 

where E + = E + iO, $i(pi) = 0i(xi)e lpi ' yi is the solution to the channel Schrodinger 
equation 

{Ho + Vf -E)$i(pi) = 0, (9) 

and Gi is the channel Green's function G\(z) = (z — H — Vf) -1 . The LSE (jSJ) is the 
integral equation of the form 

M/ + (X, Pl ) = $ 1 (X, Pl ) + 

J dX' G X (X, X', E+)[V 2 c (x! 2 ) + Vf (x' 3 ) + \gW 2 (x! 2 )}^ + (X\ Pl ) (10) 

where X = {xi,yi} and x' 2 , x' 3 are supposed to be represented through x' 1; y[ by 
standard transformations of Jacobi coordinates. The asymptotics of \l/ + (X,pi) as 
2/i — > 00 can easily be evaluated now from ffTO]) by taking the asymptotics of the Green's 
function 



Gx^X^oc-^Mxi)- ^(X'.^-eiyi). (11) 

27m J/l v 

As the result, we get ((71) with the following expression for the amplitude F 

F(pi,Pi) = f^i(p'i)l^ + ^ + W 2 |* + ( Pl )> (12) 
lira 



e WE-et y\ 
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where the matrix element stands for the integral 
($i(p'i)I^> C + V 3 C + kjW 2 |* + (pi)> = 

/ dX'$;(X',pJ)[lf (x£) + Vf K) + i^ 2 (x' 2 )]*+(X' )Pl ). (13) 

The formula (|12p suggests that the scattering amplitude can be split into the sum of 
two terms 

F = F° + igF\ (14) 

where F° is exclusively due to the Coulomb interactions between the positron and 
the hydrogen and igF 1 is due to the absorption potential only. Let us note, that the 
immediate identification of these amplitudes with pieces of ( fT2l) does not lead to the 
consistent form of the amplitudes F k since the wave-function itself may be split into 
two parts similarly to ( fl4|) . Hence, the contribution from different kinds of interactions 
cannot be separated on the basis of equation (fl2|) . The appropriate way is to rewrite 
the LSE (IE]) in the form of distorted wave representation [14] . To this end, let us recast 
(ED into 

[/ - G 1 {E + ){V 2 C + Vf p + ( Pl ) = $! + G 1 (E + )igW 2 ^ + ( Pl ). (15) 

If the energy is below the positronium formation threshold, the inversion of the operator 
from the left hand side can be performed with the help of the formulae 

[/ - G 1 {E + )(V 2 ° + Vf )]-% = ^ 0+ , (16) 

[/ - G 1 {z){V? + VfT^z) = G c (z). (17) 

As the result, the LSE (JHJ) takes the form 

# + ( Pl ) = ^ 0+ ( Pl ) + G c (E + )igW^ + ( Pl ). (18) 

The asymptotic analysis of equation ( fl~8j) will give us the representations for amplitudes 
F° and F 1 . 

The inhomogeneous term \I/ 0+ (pi) is the outgoing solution to the e + — H scattering 
problem without an absorption potential . Following ffl6|) , this function is defined by 
the solution of the LSE 

*°+( Pl ) = ^(pO + G X {E + ){V 2 C + K 3 c )*°+( Pl ). (19) 

Similarly to ([H]), the asymptotics of the solution has the form 



e VB-£i 2/1 



* 0+ (pi) oc M*i)ie lpi ' yi + ^Vyi,Pi)] (20) 

yi 

with the amplitude given by 

np'i,pi) = ^(®M)\v 2 ° +v 3 c \*° + ( Pl )). (21) 

The Green's function G (z) in (fl8|) is defined as G c (z) = {z — H Its asymptotics 
as yi — > oo reads 



_ 7/l pWE-tx yi 

G C (X,X',£+) ex * -*(X' )V /£-eiyi)- (22) 

2ixhr yi v 
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Introducing this asymptotics into (TT51) we get the explicit representation for the 
amplitude F 1 from ffT4"j) 

^(Px,Pi) = ^(^(J&m^fa)). (23) 

In these formulae \I/ 0— (pi) is the solution of the e + — H scattering problem with incoming 
boundary conditions \I/ 0— (pi) oc [\l/ 0+ (— pi)]* when only Coulomb interactions are taken 
into account in the Hamiltonian. The formulae (12TT) and fl23|) provide us with the desired 
representation for F as the sum of two amplitudes, one of which (F°) is exclusively due 
to the Coulomb interactions between the positron and the hydrogen and does not depend 
on the absorption potential and the other one (F 1 ) is due to the absorption potential 
only. 

The preceding analysis is not applicable if the energy E is higher than the 
positronium formation threshold. Indeed, whereas the formulae (JMH1) remain valid, 
the inversion of the operator / — Gi(E + )(V 2 c + V$ p) in the left hand side of (fT5|) cannot 
be performed and, as the consequence, the equations ( fT8ll23l cannot be justified. The 
formal reason is that the homogeneous equation 

X = G 1 (E + )(V 2 C + V 3 C ) X (24) 

now possesses the nontrivial solution x — ^ 02 (P2) such that 

(H c -E)V 02 (p 2 ) = (25) 

with the asymptotics \l/ 02 (p 2 ) oc 2 ( x 2)e lp2 ' y2 . The latter describes the scattering of the 
proton off the positronium ground state 02 (x 2 ) with the binding energy e 2 - 

This problem with LSE is well known in the three-body scattering theory ^5| [T6 | [T7] 
and is the manifestation of the general fact that no single LSE specifies the three-body 
scattering wave function uniquely, if the rearrangement channel is open. The resolution 
of the problem has been found by transforming LSE into the matrix equations for 
the components of the wave function. Proper arranging of the interactions between 
the equations, which guarantees the uniqueness of the solution, leads to the equations 
known as the Faddeev three-body equations [16], [TFJ [IE]. We adopt this formalism 
to our case in the next subsection where it is proven that for the energy above the 
positronium formation threshold, similar to (|T4l) . the amplitude F is given by the formula 
F = F®! + igF^ where the amplitudes F^ and F^ will be defined by formulae (15T1) and 

§3). 

2.2. Three-body scattering formalism above the positronium formation threshold 

Let us emphasize that the scattering problem with the Hamiltonian ([1]) always deals 
with the single-arrangement channel e + — H due to the presence of an absorption 
potential. That means the solution to the Schrodinger equation ([6]) has the single- 
arrangement asymptotics ([7j) irrespective that is the energy below or above the Ps- 
formation threshold. However, in order to specify the amplitudes F° and F l one needs 
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the solutions to SE with the Hamiltonian H c . For the energy above the Ps-formation 
threshold between Ps(n = 1) and H(n = 2) thresholds the Scrodinger equation 

H c m Ql = Em 0i (26) 

has two kinds of solutions, which are specified by the asymptotics 

tt^Pi) cc <Mx,)[e ip ^ + — /«], Vi -> oo, (27) 

Hi 

^ 0t (pi) ex 0fc(x fc ) /«, 7/^ > oo, k^i. (28) 

Here the momenta p^) are related to the energy as E = h 2 /2v^ k ) p?, k \ + e^) and 
indices i(k) run over the {1,2} set. The formulae (I27T [28]) reflect the fact that now the 
asymptotic form of the wave function is different in different asymptotic arrangements. 
This is exactly that property of the three-body wave-function, which cannot be recovered 
by any single LSE. 

In order to take into account the multichannel character of the scattering problem 
above the rearrangement threshold we use the formalism of Faddeev equations [18J. 
Since the original formalism is developed for the short range interaction, at the first stage 
we reformulate the three-body Hamiltonian in such a way that Coulomb interactions 
are split into long-range and short-range parts [19] 

Vf(x l ) = ^(x i ,y i ) + ^(x ij y i ). 

This splitting is made in the three-body configuration space by a smooth splitting 
function £j(xj,yj) constructed such that (^(x^yi) = 1 if Xi/xo < (1 + yi/yo) u and 
Ci(xj, yi) = if Xi/x > (1 + yi/yo) u for some x > 0, y > and < v < 1/2. With 
such a Q the short- and long-range parts of the Coulomb potentials are defined as 

V-( Xi , y< ) = Ci(x i5 y ) Vf (x,); V} = V t c - V{. 
The Hamiltonian ([1]) is then transformed into 

H = H l + V{ + V 2 S + igW 2 ; H l =T + V} + V l 2 + Vf. 

After this modification, the components of the wave function are defined by formulae 

^+ = (E + - H l )- X V^ + , (29) 
= (E + - H l )-\V 2 S + igW 2 )^ + , (30) 

where E + — E + iO. Two components are enough in our case. Indeed, the potential 
V% between the positron and the proton is repulsive and does not supports bound 
states. Hence, only two asymptotic arrangements are possible, which are covered by 
components \l/^ 2 . 

It is straightforward to see that the sum of the components recovers the wave 
function 

^+ = $+ + ^+ (31) 
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and the components obey the set of modified Faddeev equations (MFE) 

(E-H l - V{)Vt =V^+, ( 32 ) 

[E — H l — V 2 S - igW 2 )^t = (V 2 S + igW 2 )^t- (33) 



The important feature of equations ( 1321 1331 . with regard to the description of the 
annihilation, is the fact that now the two-body absorption potential igW 2 is incorporated 
into equations in two manners. It appears in the diagonal part on the left hand side of 
the equation ( 1331) and in the coupling term on the right hand side. The diagonal terms 
of the equations (|32| 1331 are responsible for the asymptotic configurations in channels 
1 and 2, so the absorption potential in the left hand side of (|33l) is natural to associate 
with the annihilation after positronium formation. The coupling term in fl32| [33]) couples 
equations in the reaction volume, therefore the absorption potential in the right hand 
side of ( |33|) is the source for the direct annihilation. These associations will be put on 
the solid ground in the next section. 

As in the case of equation (jBJ), the scattering solution to equations (13211331) is defined 
by the asymptotics as yi — > oo 

pipm 

^+( Pl ) oc <Mxi)[e lpi - yi + — F n ), (34) 

*£(pi) cx 0. (35) 
Certainly, for the amplitudes F 1X and F from (Tl2|) the equality holds true 
F n = F 

as it should be in view of ( 13~TT) . Again, as in the case of equation (JBJ), the solution 
required by (J3U [35J) is given by the integral form of MFE ([321 [33]) (IMFE) 

*f( Pl ) = $? + (pi) + (E + -H l - Vi)~ 1 W(pi), (36) 
*+( Pl ) = {E + -H l - V 2 S - igW 2 y\V 2 s + i<^ 2 )*+( Pl ). (37) 

Here $5 + stands for the solution to the channel Schrodinger equation (H l +V± — E)&i + = 
0, which explicit form is 

(H + Vi- £)$?+ = -{V 2 l + Vf )$?+. (38) 

Repeating reasoning of ([ST TTB"]) we get the asymptotics as y\ — > oo 

$? + ( Pl ) oc $!( Pl ) + & (39) 

with 

&(Pi,Pi) = ^2^i(p / i)I^ + K s C |$ 1 (Pi)). (40) 



The IMFE ( 1361 1371) are proven to have the unique solution below as well as above 
the rearrangement threshold [T8"l IT9] . We use this property to remedy the shortcoming 
of the LSE approach outlined above in specifying the structure of the amplitude F above 
the Ps-formation threshold. 





" H l + V{ 
















H l + V 2 S + igW 2 _ 
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In order to proceed, it is convenient to introduce matrix notations 

H(fc) 
V 

Thus, the MFE set takes the form 

[H{ig) - M]¥( Pl ) = -[V + i<?W]*( Pl ), (41) 

where the vector solution is constructed from wave function components as VP = 
The IMFE set p)] 1571) in matrix notations reads 

*( Pl ) = $ 0+ ( P i) + [£ + I - H(i^)]- x [V + i<?W]*( Pl ). (42) 

The vector of the inhomogeneous term is defined as <fr 0+ ( P i) = ($ ( j )+ ,0). This equation 
can be reformulated in the form of distorted wave representation by similar way we 
made for LSE (jSJ). The equation (142]) can be recast into 

{I - [E+I - H(i( 7 )]- 1 V}*(p 1 ) = d>°+( Pl ) + [E+I - H(i^)]- 1 i^W*( Pl ). 

Then, by the use of the formulae 



{I - [E+I - H(i^)]- 1 V}- 1 $ 0+ ( Pl ) = *°+( Pl ), 



(43) 
(44) 



(45) 
(46) 



{I - [*I - HO^V}- 1 ^ - H(i^)]- 1 = [zl - H(ig) - V]" 1 
we arrive at the desired distorted wave representation of the IMFE 

*( Pl ) = * 0+ ( Pl ) + [E+I - H(kj) - V]-HgW*( Pl ). 
From this equation it is seen that 

*( Pl ) = *°+( Pl ) + * 1+ ( Pl ) 

where \I/ 1+ ( P i) stands for the second term in equation fj45|) . The equation fj45|) is the 
direct analog of ( {TBI but, in contrast to that, is well defined below as well as above the 
Ps-formation threshold and therefore can be used to get the necessary representation 
for the amplitude Fn. 

The asymptotics of ^(pi) is formed from contributions of both terms in (1461) . Let 
us consider the inhomogeneous term \l> 0+ ( Pl ). The equation for this term reads 

*°+( Pl ) = <&°+( Pl ) + [E+I - H(i^)]- 1 V* 0+ ( Pl ) (47) 



or 



[H(i^ + v]* o+ ( Pl ) = M* 0+ ( P i; 



(48) 

These equations are quite similar to (1321 1331) and (1361 1371) except the coupling term V, 
which does not contain the absorption potential. The asymptotics of the solution to 
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as yi — > oo, which follows from (1471) . is similar to (1541 155^) and for the components 
of * 0+ = (^? + , ^° 2 + ) has the form 

tf^(pi) oc 0i(xi)[e lpi ' yi + — —Fix] (49) 
^ + (pi) oc 0. (50) 

As it was demonstrated above, the formula for the amplitude should be derived by 
taking asymptotics of the right hand side of ( 14"T1) as y\ — > oo. A minor difference from 
what we have demonstrated above is that the nontrivial contribution into the amplitude 
comes from not only the Green's function source term but also from the driving term 
due to (139|) . So that, the amplitude F^ is given by 

F °i (p, i' pi) = ^[(^i(p;)i^+^i^r(pi))+<^r(p;)in s i^2 + (pi)>]-(5i) 

Now we consider the ^ ,1+ (pi) term. Its asymptotics as yi — > oo is defined by the Green's 
function G(z) = [zl — H(ig) — V] -1 . This function is a genuine three-body quantity 
and the asymptotics of its matrix elements as yi — > oo may be written in the form 



p l\/E-eiyi 

G ij (X,X',E+) oc — ^(x 4 ) T<j-*(X!,y/E-e i y i ), (52) 

Z7T/T yi 

where T°~ is the eigenfunction of the adjoint to the operator from the left hand side of 
( 1481) . This asymptotics has different character for i = 1 and i = 2. In the first case G\j 
does not vanish with yi large since E — ei is real nonnegative. If i = 2, the Positronium 
binding energy e 2 becomes complex when the absorption potential is introduced into 
the e + e~ Hamiltonian. That makes the relative momentum p 2 = y/E — e 2 complex, 
i.e. p 2 = p 2 + ip\. So, the asymptotics of G 2 j vanish exponentially. Thus, only \V{ + 
component has the nontrivial asymptotics 

gipis/i 

ttl+CpiWiCxO— —igFli (53) 
with the amplitude Fn given by 

^i(pi,Pi) = ^f<T^(pi)|W 2 |*f(pi)>. (54) 
Here T 2 ~ is the second component of the solution to the adjoint equation to (j48H 

[H(-i^) + V T - J E;i]Y -(p 1 ) = (55) 
being defined by the integral form 

Y°-( Pl ) = *°-( Pl ) + [E-I - H(-i^)]- 1 V T T°-(p 1 ). (56) 

In this equation E~ = E — iO, V T means transposed matrix and <fr 0- = ($? _ , 0) where 
is the solution to (1381) with incoming boundary conditions. The solution to (1551) 
in the case of g = takes the very simple form, i.e. = T 2 = \V° where \V° is the 
three-body wave-function for the pure Coulomb problem. So, \V° obeys the Schrodinger 
equation (|26|) and can be constructed from wave-function Faddeev components (|29| [30j) 
at g = as \V° = \Vi° + \V 2 °. The detailed information on this and other features of the 
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matrix equations as MFE and IMFE and matrix Green's functions, which is necessary 
for evaluations made above, can be found in Ref. [20J. 

The formulae ( |5Tl 15^1) determine the ingredients of the amplitude Fu uniformly 
above as well as below the rearrangement threshold of the positronium formation. It is 
possible to show by preforming backwards transformations that below the Ps-formation 
threshold, where the representations (TH] I2TI 123]) are valid, the equality 

F* 1 + igFl 1 = F° + igF 1 = F (57) 

holds true. It is important to note, that from the analysis of the equation (HTj) it follows 

Ff 1 = F + O(g). (58) 

This means that the respective terms in the left and right hand sides of the first equation 
in the chain ( 1571) are not identical. In fact the amplitude F^ takes into account not 
only the Coulomb interaction between the positron and the hydrogen, as F° does, 
but also the possibility of e + e~ annihilation after the positronium formation. That 
is due to the presence of the absorption potential in the diagonal part of the equations 
( |48|) . This means that below the positronium formation threshold the annihilation 
after virtual formation of the positronium is incorporated into the F{\ amplitude and 
consequently F^ is the pure direct annihilation amplitude. It is worth mentioning again, 
that the formulae ( 1511 15^1) determine the amplitudes uniformly below as well as above 
the positronium formation threshold, whereas it is not true for representations (12TI 123]) . 
They are valid only below the positronium formation threshold. 



3. Optical theorem in the presence of absorption and annihilation cross 
section 

3.1. Optical theorem 

The standard optical theorem for Hermitian Hamiltonians is nothing but the 
manifestation of the flux conservation, what is equivalent to the unitarity of the S- 
matrix. The absorption potential breaks the Hermiticity and the scattering is not 
unitary. The lack of unitarity is the measure of how much of the flux is absorbed and 
in the case of the annihilation is the way to determine the annihilation cross section. 
There is extensive literature on the optical theorem but [H] is the most suitable for our 
purpose. Following this approach by multiplying the SE (Q by complex conjugate wave 
function and subtracting the complex conjugated SE multiplied by we arrive 
at the equality 

^ + *H ^ + - ^ + H ^ + * = -2igW 2 \^ + \ 2 . (59) 

Integrating over the domain = {yi < R}, using the Green's formula and taking the 
limit as R —>■ oo, we get the following result 

- 2ig / d^dy, W 2 \^ + \ 2 = ^- lim / R 2 dy! f dx a {^+, ^+*}, (60) 
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J d^d yi W 2 \^\ 2 = — 3mF( Pl , Pl ) -Jdy |F( Pl y, Pl )| 2 . (61) 



which represents the balance of the flux. Here the Wronskian ^ + *} = \E r+ <9 2/1 \E r+ * — 
^ + *d yi ^> + has to be taken at the condition y\ = R. Then, using the asymptotic form 
of given in (CO), the normalization of the ground-state wave functions 0i(xi) and the 
weak asymptotics of the plain wave (see for example [2T]) 

e W oc — \-5(p + y)e Apy + 5(p - y)e ipy 
ipy ' 

we finally arrive at the optical theorem in the presence of absorption 

n pi J pi 
The positive quantity in the left hand side of (IBTj) is the absorption cross section which 
determines in our case the overall e + e~ annihilation cross section due to the direct 
process as well as due to the annihilation after (virtual if E < e 2 or actual if E > e 2 ) 
positronium formation 

o* = ^t^> f d x ldyi W 2 \V+\ 2 . (62) 
n pi J 

So that, the overall annihilation cross section can be computed either by the integral 
( I621) or by the expression in the right hand side of ( l6Ti) . if the total amplitude F is in 
possession. 

In order to go beyond the standard formulation (l6Tl 1621) one needs to use the 
detailed structure of the amplitude F. In our case it is the representation (1571) 

F = F° + igFl, 

which leads to the following form of (16T1) 

a a = a a 2 + o\ (63) 

a a 2 = —%mF? 1 ( Pl , Pl ) - / dy\F^( Pl y,Pi)\ 2 (64) 
Pi J 

Att 

0-1=— (-^SteFiGh.pi)- 
Pi 

2{-g) J dyZmF^p^pjF^ipiy^-g 2 ! dy {F^foy, Pl )| 2 . (65) 

The quantities af and have meaning of annihilation cross sections for the direct 
process and the process of the annihilation after the positronium formation. To make 
this statement sound, let us show that a 2 is the cross section of the after the positronium 
formation annihilation. Thus, the remaining part of a a , what is <r", should be interpreted 
as the direct annihilation cross section. 

It is apparent, that the equation fl64l) is the optical theorem formulated for the 
equation (jUJ). Indeed, multiplying the equation (jUj) by T +0 from the left and 
subtracting the equation for Y +0 

[H(i2)+V T -£I]T 0+ ( Pl ) = ( 66 ) 

multiplied by \I/ 0+ from the right and making obvious cancelations we get 
( Y 0+ , H W 0+ ) - (H ir 0+ ,V 0+ ) = 2h;(Y 0+ ,M 0+ ). 
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a 



2 _ fe 2„ 



I dx 1 rfy 1 T° + *(p 1 )^ 2 vI/° + (p 1 ). (6? 



Here D is a diagonal matrix diag{0, W 2 } and (., .) is the scalar product in the two 
dimensional complex space C 2 of wave-function components. Repeating argumentations 
which led us to the formula ( 16T1) . we arrive at the equality 

2 -^- / dx^rrMw^Pi) = 
a pi j 

— Qtoi^CPLPO - / rfy I^C^y, Pi)| 2 - (67) 

This is the optical theorem for equation fHHj) and the annihilation cross section a% can 
be expressed now in terms of the left hand side as 
2vj(-g) 

h 2 Pl 

To elucidate the further meaning of the cross section o~% it is instructive to consider 
the limiting case as g — > 0. In the limit no absorption potential is present, scattering 
becomes unitary and amplitude coincides with / n from (127|) . The standard unitary 
variant of the optical theorem for the amplitude / n has the form 

^Ota/u^pO- / rfy|/n(Piy,Pi)| 2 - / rfy|/ 2 i(Piy,Pi)| 2 = o (69) 
Pi J J 

where fu is the elastic e + — H amplitude and f'21 is the rearrangement Ps— p amplitude. 
Therefore, the right hand side of fl67|) has the limit 

An, 
Pi 

The quantity o~ii is nothing but the positronium formation cross section. At the 
same time the formula fl68|) leads to an uncertainty in the limit as g — > when the 
diverging integral is multiplied by the vanishing factor g. Actually, the resolution of 
this uncertainty is made by flTUj) and gives the relation 

a« = a 2 i + 0(g), (71) 

which clearly shows that is the cross section of the annihilation after the positronium 
formation. 

The main result of this subsection is the representation of the annihilation cross 
section as the sum of two terms 

a = a x + cr 2 , 

where is shown to represent the cross section of the annihilation after the positronium 
formation, which is given by the formulae (1M1) or (16"S|) . Therefore a" is the direct 
annihilation cross section, given by (|65l) . 

The analysis made above shows that the definitions of cross sections in terms of 
scattering amplitudes fl64t 165]) are uniform and are valid below as well as above the 
rearrangement threshold of the positronium formation. Moreover, by construction, the 
integrals ( 15T1 154}) involved in the definitions of the amplitudes F X1 , F X1 have the finite 
limit as g — * 0. This property will be used below for perturbative calculations of the 
amplitudes and cross sections. At the same time the integrals in representations (1521 1681 ) 



— Qtai^pLpO- / dylF^y.pOl 2 ^ f dy\f 21 ( Pl y, Pl )\ 2 = a 21 .(70) 
pi J J 
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are divergent above the positronium formation threshold if the limit g — > is taken. 
That makes these formulae not suitable for perturbative methods above the positronium 
formation threshold. 

3.2. Absorbing annihilation potential 

In this subsection we fix the coordinate form of the absorption potential by comparing 
the definition (1621) for o a with QED formula for the 27 singlet e + e~ annihilation |jj [21 E] 

a a = nr 2 (c/v)Z eil . (72) 

In this formula r$ is the classical electron radius, c is the speed of light and v is the 
incident velocity of the positron. The effective number of electrons Z g g participating in 
annihilation is given by the integral 

Z eS = J l^° + (xi,yi)| 2 5(x 2 ). (73) 

Here \l/ 0+ is the \l/ 01 solution of the e + — H scattering problem ( |26l [271 BHD w hen the 
absorption potential is not taken into account. The integral in (1731) is well defined below 
the Ps-formation threshold. In this case (1721 1731) can be considered as the first order 
perturbation approximation to fl62l) . since below the rearrangement threshold ~ \l/ 0+ 
is the well defined first order perturbation solution (Born) to ( |T8|) . This observation was 
used in [I] to determine the absorption potential for 27 singlet e + e~ annihilation as 

igW 2 {* 2 ) = iao(x 2 ) (74) 
e 2 

g = 2na 3 . 

a 

Here ao is the Bohr radius and a is the fine structure constant. This potential was used in 
[I] to calculate the direct annihilation cross section by solving the Lippmann-Schwinger 
equation for T-matrix below the Ps-formation threshold. 

The formula (!73|) cannot be extended for calculations above the Ps-formation 
threshold since the integral diverges. One of ways to go over the Ps-formation 
threshold is the use of (1521 with absorption potential (1741) incorporated into the 
Schrodinger equation. This was done in papers [5J [7] where the overall annihilation 
cross section a a was computed below as well as above the Ps-formation threshold 
within the hypersperical close coupling technique for the time independent three-body 
Schrodinger equation. Paper [8] represents the solution of the time- dependent three- 
body Schrodinger equation with the absorption potential (I74|) below as well as above the 
Ps-formation threshold. All these papers dealt with numerical solutions of respective 
equations and the delta-functional singularity of the potential ( 174|) was treated by a 
certain numerical approximation. Nevertheless, the analytical status of the potential 
( 1741) is not satisfactory. The delta-functional singularity is too strong and makes the 
Hamiltonian not well defined. This issue was not addressed in aforementioned papers 
and we give a portion of the necessary analysis in this subsection. 

% The spin-averaging factor 1/4 for singlet 27 annihilation is implied implicitly. 
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It is well known already since the papers by Fermi [5] and then Breit [TD] that 
the three-dimensional delta-function potential can be incorporated into the Schrodinger 
equation only perturbatively. One of the approaches to go beyond the perturbative 
treatment is the use of a zero-range potential [TOl [TT] . There are two common ways to 
introduce the zero-range potential. One is imposing boundary conditions for the wave- 
function at the singularity point. The other one is introducing into the Hamiltonian 
an additional term, which enforces the wave-function to fulfill the boundary conditions. 
This term can conveniently be represented in the compact form by the quasi-potential 
[22| [23] . We choose the second option. 

The singularity caused by the zero-range potential in the case of the electron- 
positron interaction is located at the same point as the Coulomb singularity —e 2 /x 2 . 
The latter leads to the modification of the standard zero-range potential and of the 
respective quasi-potential. The resulting definition for the coordinate part W 2 of the 
absorption potential is 

W 2 (x) = g(x) l ~ n f ± X - (75) 

1 + n 2 x log x dx 1 + n 2 x log x 

where n 2 = — 2/i 2 e 2 and /i 2 is the e + e~ reduced mass. The detailed derivation of (1751) 
involves a substantial portion of mathematics and will be published elsewhere. Some 
basic theorems, which define the general properties of the zero-range potential with the 
Coulomb modification, can be found in [2lj . 

It is straightforward to see that in the limit n 2 — * the quasi-potential (1751) takes 
the standard form 

W 2 (x) ->$(x)-^-x. (76) 
dx 

The quasi-potential igiy 2 (x 2 ) enforces the following asymptotics for the wave function 
as x 2 —* 

^(x 2 ,y 2 ) oc + n 2 \ogx 2 }+b(y 2 ) + 0(x 2 \ogx 2 ), (77) 

4:71 X 2 

where a/b = —ig2fx 2 /h 2 . This asymptotics, as usually, determines the appropriate 
boundary conditions, which we do not write down here explicitly. It can be shown that 
the action of the quasi-potential on the function with such an asymptotics is given by 
the formula 

^ 2 (x 2 )^(x 2 ,y 2 ) = <J(x 2 )6(ya). (78) 

The latter means that the action of the quasi-potential on a function x(x 2 ,y 2 ), which 
is regular in the point x 2 = 0, is equivalent to the delta-function 

W 2 (x 2 )x(x 2 ,y 2 ) = <5(x 2 ) X (0,y 2 ). (79) 

This formula shows that if a matrix element of the quasi-potential is calculated between 
the functions \i <^> which are regular at x 2 = 0, then the quasi-potential is equivalent 
to the delta-function, i.e. 

(x\W 2 ( X2 )\lu) = (x\5(x 2 )\cu). (80) 
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0.5 
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0.061 
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0.6 


0.0032 


0.003 


0.0039 


0.003 


0.003 


0.7 


-0.0502 


-0.05 


-0.051 


0.0053 





Table 1. L = phase shift So for e + — H elastic scattering. Momenta pi are given in 
units of (Iq 1 . 

This statement justifies the use of the delta-function as the coordinate part of the 
absorption potential in [51 [71 [8], since the basis functions of the approaches used to 
compute the matrix elements of the absorption potential are smooth. Nevertheless, 
any basis of smooth functions cannot reproduce the singularity in f!77|) by a finite 
number of terms, what always happens along the numerical solution. Hence, such a 
treatment of the zero-range potential in [2J UJ [8] is approximative but, in view of the 
fact a(y2) = —ig2^2o{y2)/h 2 with \g\ <^ 1 for the e + e~ annihilation, the approximation 
is reasonable. 

4. Calculation of annihilation cross section 

The formalism developed above was applied for calculations of the annihilation in e + — H 
collision below the Ps-formation threshold and above the threshold in the Ore gap. The 
latter is defined as the interval of the energy between the Ps(ls) and H(n = 2) thresholds. 
All calculations were made on the basis of an extension of the multichannel numerical 
algorithm for Faddeev equations, described in details in [23] , for the case of annihilation. 
The algorithm uses the bipolar harmonic expansion to represent the angular dependence 
of the wave-function components 

^ + (x,y.) = E ^ h(z< ' y,) y^,yO, (81) 

Lhi 2 Xr >-y % 

This expansion reduces the MFE to a set of coupled equations for radial components 
i>i x i 2 {xi, yi) which are then approximated by the quintic-spline expansion and solved by 
the orthogonal collocation procedure. The maximum values of h,fa used in (18TT) range 
from 12 to 15. 

In order to test the numerical approach, the genuine Coulomb problem for e + — H 
scattering was solved on the platform of MFE (|32| [33]) with g = 0. The results of 
calculations for s-wave e + — H phase shift So are given in Table [U together with data of 
other authors. As one can see, the agreement of our calculations with previous results 
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Figure 1. L = phase-shift for e + — H clastic scattering. 



is quite good. Figure [T] provides an alternative representation of 5 as the function 
of the energy to demonstrate the regular character of the calculated phase shift 5o- 
The calculations of the direct annihilation cross section er" were performed with the 
help of the representation fl65l) . The amplitudes and F^ were computed from the 
solutions of the respective Faddeev equations. The effective number for the direct 
annihilation cross section was obtained using the standard expression (IT2]) 

cr? = 7rr 2 (c/v)Zl s . (82) 

As it was mentioned in subsection 13. 1[ in the calculations of the amplitudes F^ and F^ 
we have systematically approximated the Faddeev components, involved in the matrix- 
elements for the amplitudes, by the solutions \E , ° 0± of the MFE with g = 0. This 
statement can easily be justified by the iterative solution of the IMFE ( 136| 1371) for 
| g | <C 1. The iterative solutions is well defined thanks to the fact that the matrix 
kernel of the Faddeev equation is proven to be compact. This is another advantage 
of the Faddeev three-body equations. As the result, the expression for the absorption 



amplitude from ([5 
*ii(Pi,P) 



can be simplified as 



([*r(Pi)+*r(Pi)]iw 2 i^ u+ (pi)) 



2nh z 



00+/ 



(83) 



which was actually used for calculations of this paper. 

Since a number of data from other authors for direct annihilation cross section below 
the Ps-formation threshold is available, In Table [2] we display our calculated phase-shift 
So and Z^g together with results of other authors obtained with the standard formula 
(172]) for one of the typical value of the relative momentum pi = 0.4 [1/ao] • The agreement 
is very good for calculations made by quite different approaches. In Table [3] we collect 
the results of existing calculations of Z g g below the Ps-formation threshold to compare 
with our results. One cannot expect the complete agreement since our definition of a\ 
concerns the direct process only, whereas the standard definition below the Ps-formation 
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Table 2. L = phase shift Sq for e + — H elastic scattering and effective number Z & q 
for the relative momentum p\ — 0.4 [a^" 1 ]. 
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2.4625 


2.29 


2.3 
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0.7 


2.2529 


2.02 




1.950 



Table 3. L = effective number Z q q for e + — H annihilation. The relative momenta 
Pi are given in units of a^ 1 . 

threshold deals with the overall annihilation cross section a a . Nevertheless, Table [3] 
shows that the difference is not so dramatic. 

The extensive calculations of the direct annihilation cross section were made in 
the interval of the energy between H(n = 1) and H(n = 2) thresholds on the basis of the 
formula (|65j) . Figure [2] shows the s-wave effective number derived from af, with the 
formula (18"2"|) . and the results of calculations from [8j. Although, quite different methods 
(time-independent Faddeev equations and time-dependent wave-packet approach to the 
Schrodinger equation) of calculations for the direct annihilation cross sections have been 
used, the agreement between data is fairly good. In fact, the definition of the direct cross 
section cr" in our formalism as the remainder of the overall annihilation cross section o~ a 
after subtraction of the annihilation cross section after the Ps-formation <jf 

a{ = a a - a a 2 , (84) 

(see equation ( |63l) ) is quite equivalent to the time-dependent definition by formula (7) of 
ref. [S]. Nevertheless, as the data of Tables EE \5\ shows, the exact use of the formula (1841) 
would be very unpractical. The formation cross section cr 2 i, which is the leading term of 
o"2 due to (ED) , rapidly increases above the threshold and is several order of magnitude 
bigger than cr®. The analytic separation of the overall annihilation cross section into 
the formation and direct parts made in Subsection 13.11 and given by formulae ( 16311651) in 
terms of amplitudes and is therefore of the great practical importance. 
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Figure 2. The effective number for the L = direct annihilation. The diamonds 
are the calculations of the present paper, the open squares are the data from [5], the 
triangle marks the positronium formation threshold. 

In the following Tables H] and we present the results of our calculations for the 
direct annihilation cross section together with the cross section of the positronium 
formation above the Ps-formation threshold. They are given for slightly different 
values of momenta p\ for L = and L = 1 in order to emphasize the most characteristic 
behavior of cross sections near the respective Feshbach resonances. Besides the expected 
difference in several order of magnitude between direct annihilation cross section cr" and 
the positronium formation cross section cr 2 i, the strong correlation between these cross 
sections in the region of the sharp increase of cr" across 0.86124 [a^ 1 ] for L = and 
0.86318 [qq 1 ] for L = 1 is clearly seen. The graphical representation of that sharp 
increase of cross sections is given on Figures [3] and H] which display s- and p-wave Z^q. 

This resonant feature was also observed in the spatial behavior of the solution to the 
Faddeev equation. The first component tyf , which is related to the e + — H channel, 
always exhibits the characteristic resonant bump near resonant energy. These spatial 
resonant structures can be extensive depending on the lifetime of the resonances. In 
[T3] , reasonable estimation of the energy width of a number of Feshbach resonances were 
made using the uncertainty principle applied to the graphical representation of the wave 
function components. 

5. Conclusion 

In present paper we have developed the multichannel time-independent formalism, which 
is capable to describe the scattering and annihilation processes in the positron-hydrogen 
collision above the rearrangement threshold. The expression derived in the paper for the 
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Table 4. L = effective number •^ 1 g-, the direct annihilation cross section af, and 
the positronium formation cross section ai\. The cross sections are given in units of 
irciQ and momenta in units of Oq 1 . The abbreviation [-n] is used for 10 _n . 
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Table 5. L = 1 effective number Z^, the direct annihilation cross section af, and 
the positronium formation cross section a 2 i- The cross sections are given in units of 
7rao and momenta in units of Oq 1 . The abbreviation [-n] is used for 10 _n . 
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Figure 3. L = effective number in the resonant region. The momenta p are 
given in a^ 1 units. 
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Figure 4. L — 1 effective number in the resonant region. The momenta p are 
given in a^ 1 units. 
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direct annihilation cross section in terms of amplitudes is proven to extend the standard 
formula to the energy region above the positronium formation threshold. Below the 
threshold our cross section and the cross section calculated from the standard theory are 
in good agreement provided the close vicinity of the positronium-formation threshold is 
not considered where the standard cross section becomes infinite. The direct annihilation 
cross section defined by the formula fl65l) does not exhibit any singular behavior at the 
Ps-formation threshold. It is in good agreement with the nonsingular direct annihilation 
cross section computed from the time-dependent solution of the three-body Schrodinger 
equation for e + — H system [8]. 

The formalism of this paper can readily be extended beyond the Ore gap. In this 
case the multichannel optical theorem, which generalizes fl61|) . plays the key role in 
the determination of cross sections. Preliminary calculations indicated the much larger 
enhancement of the direct annihilation cross section near the Feshbach resonances in 
the eight open channel region above the Ps(n = 2) threshold. Since there are numerous 
resonances beyond the Ore gap, they evidently made a significant contribution to the 
overall annihilation peak around 12 e.V. displayed in Figure 2 of [8]. We hope that 
our individual annihilation resonance structure will provide a new experimental tool to 
study sharp resonances. Although at the present time, it is not possible to conduct such 
an experiment for e + — H system, experiments for positron scattering on large molecules 
have been done for many systems [34J. 

Our approach, in perspective, opens the way to considering more complicated 
systems with more than three particles, as for example e + e + e~e~. Suitable formalism 
for multichannel scattering [35], which is the generalization of the Faddeev equations 
for the four-particle systems, will be helpful as an important theoretical step towards 
the experimental verification and utilization of the rich positron annihilation physics. 
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